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Abstract 

Nonlinear reaction-diffusion systems admit a wide variety of spatiotemporal patterns or 
structures. In this lecture, we point out that there is certain advantage in studying discrete 
arrays, namely cellular neural/nonlinear networks (CNNs), over continuous systems. Then, 
to illustrate these ideas, the dynamics of diffusively coupled one and two dimensional cellular 
nonlinear networks (CNNs), involving Murali-Lakshmanan-Chua circuit as the basic element, 
is considered. Propagation failure in the case of uniform diffusion and propagation blocking in 
the case of defects are pointed out. The mechanism behind these phenomena in terms of loss 
of stability is explained. Various spatiotemporal patterns arising from diffusion driven insta- 
bility such as hexagons, rhombous and rolls are considered when external forces are absent. 
Existence of penta-hepta defects and removal of them due to external forcing is discussed. 
The transition from hexagonal to roll structure and breathing oscillations in the presence of 
external forcing is also demonstrated. Further spatiotemporal chaos, synchronization and size 
instability in the coupled chaotic systems are elucidated. 
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1 Introduction 

Patterns abound in nature. If we look around us, we can observe myriads of interesting patterns 
ranging from uniform to very complex varieties. They occur in varied phenomena encompassing 
physics, chemistry, biology, social dynamics, economics and so on. Essentially they are distinct 
structures on a space-time scale, which arise as a collective and cooperative phenomena due to 
the underlying large number of constituent systems. The latter could be aggregates of particles, 
atoms, molecules, cells, circuits, defects, dislocations and so on. When these aggregates can move 
and/or interact, they give rise to the various patterns. A small select set of patterns is shown in 
Fig. |l|. The patterns tell us much about the dynamics of the macroscopic as well as to some extent 
the microscopic behaviours of the underlying systems. Naturally when the interactions among the 
constituents are nonlinear, one might expect novel and unexpected patterns. 

Patterns could be simple and complex. They could be stationary (eg. still image) or changing 
with time (eg. recurrent image). They could also tend towards a goal or target asymptotically. 
Homogeneous or uniform patterns, though trivial, are important basic structures. One can have 
often travelling wave patterns, especially in dispersive systems. Under suitable nonlinear forces, 
dispersive systems can even admit a different kind of waves namely solitons, which are novel type 
of localized spatiotemporal patterns, retaining their identity for ever. Perturbations of them can 
also give rise to further interesting structures. However, even more novel structures which can 
mimic naturally occurring patterns arise when one considers nonlinear diffusive (especially the so 
called reaction- diffusion) systems. When large aggregates of microstructures consisting of atoms, 
molecules, defects, dislocations etc are able to move and interact, the evolution of the concentration 
of the species can be shown to obey nonlinear diffusive equations of reaction type. They may be 
deduced from the underlying mass, energy, momentum, etc, balance equations. The general form 
of the nonlinear reaction-diffusion equation can be given by 

^^=V.{D^c) + f{c,r,t). (1) 

Here c represents the population or concentration density of the species and D and / are, in general, 
nonlinear functions of c representing the diffusivity and the reaction kinetics, respectively. 
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Figure 1: A select band of patterns in natural systems: (a) Zebra (b) Panther (c) snow leopard 
(d) peacock (e) giraffe (f) pilar coral and (e) daisy coral 



For example, one can consider the population density of a particular species and then / in 
eq. (0) will represent the birth and death processes. In the case of logistic population growth 
/ ac(l — 6c), where a is the linear reproduction rate and S is the inhibition rate. If D is a 
constant, then, the above equation (Q) will correspond to the celebrated Fisher equation, 

dc 

— = DV^c + ac(l-Sc). (2) 
ot 

Originally Fisher proposed the one dimensional version of the equation in 1937 for the spread of 
an advantageous gene in a population[0. 

1.1 Reaction-diffusion systems 

Let us consider the situation in which two or more number of species such as different chemicals, 
bacteria, populations and so on interact and evolve. In such a situation, one may study the un- 
derlying dynamics by considering the appropriate reaction-diffusion equations. Examples include 
the following systems. 

1.1.1 Belousov-Zhabotinsky reaction (Oregonator model) 

This is a simple model introduced by Fields, Koros and Noyes of University of Oregon, U.S.A. 
in 1972 to explain the various features of the Belousov-Zhabotinsky reaction (for details see for 
example [^). It can be expressed in a particular form as 



bu2iui~a) 

Ul(l - Ul) 



Ul 

U2t = D2V^U2 +Ui- U2. (3) 

Here ui represents the concentration of the autocatalytic species HBr02 and U2 is the concentra- 
tion of the transition ion catalyst in the oxidized state Ce^^ or Fe'^+. -q, a and b are parameters. 
In eq. (^) the suffix t represents partial derivative. 
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1.1.2 Brusselator model 

This is one of the often studied model for the formation of chemical patterns introduced originally 
by Lefever, Nicolis and Prigogine (see ||]). It is based on the following chemical reaction: 

A — > X, 

B + X — ^ Y, 

2X + Y — > 3X, 

X ^ E, (4) 

where the concentration of the species A, B and E are maintained constant and are thus real 
control parameters of the system. After appropriate scaling, the evolution of the active species X 
and Y can be described by the following set of equations 

dtX = A-{B + 1)X + X^Y + DxV^X, 

dtY = BX - X^Y + DyV^Y. (5) 
Here Dx and Dy are diffusion coefficients. 

1.1.3 Lotka-Volterra predator-prey model 

A model for interacting populations of two species is one in which the population of the prey is 
dependent on the predator and vice versa. Such systems can be represented by the following set 
of equations [||: 

dSi d'^Si L c c 

^ ^ D-2^-^--^S, + b,S,S,. (6) 

Here 5*1 and 5*2 are the population densities of prey and predator, Di and D2 are the diffusivities 
of the two populations, respectively, ai, 02 are the linear ratio of birth and death for the individual 
species; 61 and &2 are the linear decay and growth factors due to interaction. 

1.1.4 Gierer-Meinhardt model for biological pattern formation 

The possible interaction of an activator and a rapidly diffusing inhibitor can be modeled by the 
following set of equations Q: 

da a? 

— = Da/\a + Pa— J-Mafl + CTa, 

— = DhAh + pha"^ - ^ihh + ah- (7) 

Here fJ-h are removal rates, pa, Ph are the cross reaction coefficients and (Ja, ah are the basic 
products of the activator and inhibitor, respectively. Kq corresponds to the saturation constant. 

1.1.5 FitzHugh-Nagumo nerve conduction model 

A well known model for impulse propagation along the neuronal axons of living organisms is the 
FitzHugh-Nagumo model represented by the following set of equations 

Vt = V,, + V- — -R + I{x,t), 

= c{V + a-bR). (8) 

Here V represents the action potential and R corresponds to the lumped refractory variable. I{x, t) 
is the external injected current, a, b are positive constants and c corresponds to the temperature 
factor. 
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2 Spatiotemporal patterns in reaction-diffusion systems 



In the above, we have mentioned a few reaction-diffusion systems that arise under different phys- 
ical, chemical and biological contexts. Next it is of considerable interest to look into the distinct 
space-time structures or spatiotemporal patterns which are admitted by these systems. Specifically 
we may mention the following structures: 

• uniform or homogeneous steady states (trivial) 

• autowaves including travelling waves 

• spiral waves 

• Turing patterns (rolls, stripes, hexagons, rhombs, etc.) 

• localized structures 

• spatiotemporal chaos 

and so on. Among these patterns some of them (eg. Turing patterns) may also be identified 
as dissipative structures as they emerge spontaneously from homogeneous equilibrium states and 
correspond to systems driven away from thermodynamic equilibrium. In the following we will 
briefly discuss the various space-time structures. 

2.1 Homogeneous patterns 

The trivial, but important, class of patterns is the homogeneous or uniform steady states. These 
are the equilibrium solutions of the governing equations. Interesting dynamical features such as the 
formation of various nonhomogeneous spatiotemporal patterns will arise when these homogeneous 
states lose their stability. 

2.2 Autowaves 

Transport processes in physical, chemical and physiological systems have been often associated 
with special types of waves, namely travelling waves. In general, these types of waves are termed 
as autowaves or autonomous waves^. These autowaves can propagate in an active excitable 
medium at the expense of energy stored in the medium even in the absence of external driving 
forces. One has to note that dispersive systems often admit different kinds of propagating waves, 
called classical waves, such as sinusoidal travelling waves, wave packets, solitary waves, solitons 
and so on. However, all these dispersive waves carry energy and information and do not consume 
any energy associated with the medium. The medium in these cases is said to be passive. On 
the other hand, autowaves propagate in nonlinear diffusive as well as dissipative systems in a self 
sustained manner by inducing a local release of the stored energy and use it to trigger the same 
process in the adjacent regions. The medium here is termed as active medium. Typical examples 
include the waves of combustion (eg. flame propagating along the cord of a cracker), waves of 
phase transition, concentration waves in chemical reactions, nerve impulse propagation in neuronal 
axons, excitation waves in cardiac tissues, epidemic waves in ecological populations and so on. 

Few important classes of autowaves which arise in many physical, chemical and biological 
reaction-diffusion systems are the travelling wavefronts, travelling pulses, travelling wave trains 
and so on. An example for the travelling wavefront is the solution of the Fisher equation (|^), in 
one spatial dimension, of the form 

u — u(x ~ ct) = u{C), C, — X — ct. 

Here c is the wave speed. Travelling pulses arise for example in the FitzHugh-Nagumo nerve 
conduction model (H) . Typical example for the travelling wavetrain is the wave solution of A - w 
systems js). 
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2.3 Spiral waves 



Another interesting class of waves which are very common in a variety of natural systems corre- 
sponds to the spiral waves and scroll waves that arise in more than one spatial dimensions. The 
spiral waves which arise in two spatial dimensions can be realized by considering a single wave 
propagating around a circular obstacle The wave repeatedly travels along the same path at a 
frequency given by the wave velocity divided by the circumference of the obstacle. If the radius 
is gradually decreased the frequency of the obstacle increases until the wavefront of a new wave 
catches up with the tail of the previous wave. At this point the rate at which the wave travels 
around the obstacle cannot increase further because the new wave cannot re-excite regions that 
are still recovering from the previous wave. 

If the radius is made even smaller, the wave is forced to adopt a spiral shape that continues 
to rotate around a central core (see Fig. ||(c)). The spiral cannot enter the central core because 
this region is still in the so called refractory (or recovering) state and cannot be re-excited. An 




Figure 2: (a) Hexagonal pattern in Rayleigh-Benard system (b) Localized propagating wave in 
binary fluid convection (c) Sprial wave in the Belousov-Zhabotinsky reaction and (d) Localized 
standing wave in vertically vibrated layer of bronze balls 

example for the spiral wave is the case of cardiac arrhythmia in the heart during an irregularity 
in the heartbeat. The phenomenon is known to be modeled by Hodgkin-Huxley model. Another 
example is the spiral patterns in the Belousov-Zhabotinsky (BZ) reaction where scroll waves have 
also been observed. Recent studies show the generation of scroll waves in photosensitive excitable 
media by perturbing travelling waves to their direction of propagation pO| 

2.4 Turing patterns 

Yet another important kind of patterns which arises in reaction-diffusion systems is the Turing 
pattern. A. M. Turing in 1952 suggested that, under certain conditions, chemicals can react and 
diffuse in such a way to produce steady state heterogeneous patterns of chemicals He had 
also proposed a model for the chemical basis of morphogenesis (which represents the development 
of structure during the growth of an organism) . Turing patterns arise in many reaction-diffusion 
systems when a homogeneous steady state which is stable due to small spatial perturbations in 
the absence of diffusion becomes unstable in the presence of diffusion combined with the condition 
that one of the species or chemicals should diffuse faster than the other. Typical examples include, 
hexagonal patterns in Rayleigh-Benard system, the formation of spots in leopard, stripes in zebra, 
fingerprint patterns and so on. In order to understand the formation of Turing patterns, let 
us consider the model for the two chemical species represented by the following set of rescaled 
reaction-diffusion equations: 

— = 7/(w,w) + \/ u, 

— = jg{u,v) + dV^v. (9) 
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Here u and v are the rescaled concentration of the two chemical species, / and g are, in general, 
nonlinear functions of u and v representing the kinetics. Suppose (uo,wo) is the homogeneous 
steady state, which is stable when the diffusion is absent and becomes unstable in the presence of 
diffusion. Then one can write the condition for the diffusion driven instability asQ 

fu+9v < 0, 

fugv - fvQu > 0, 

dfu + gv > 0, 

{dfu + gvf -^d{Ugv- fvgu) > 0, (10) 

where /„, /„, gu and gy are the derivatives of the functions / and g evaluated at the steady 
state (uo. Wo)- These conditions can be further used for a selection of the admissible wavenumbers 
associated with the instability, leading to spontaneous formation of nontrivial spatial patterns. 

2.5 Localized structure 

Recently it has been found that certain dissipative nonlinear systems when driven by external 
forces exhibit a novel class of localized structures. Such structures have been observed in several 
experimental and theoretical models ||l2|, |l^, 0. For example, vertically vibrarted granuler layer 
shows a kind of localized oscillations called oscillons. These oscillons are highly localized particle- 
like excitations (see Fig. ||(d)) of the granular layer which oscillate at half the driving frequency. 
Other examples includes localized oscillations in fluid systems breathing solutions of Ginzburg- 
Landau equationsjlj] and so on. Essentially these localized structures arise due to a tendency of 
certain nonlinear sytems which localize dissipation in the presence of external forces. 

2.6 Spatiotemporal chaos 

The study of spatiotemporal or extensive chaos has also been receiving considerable interest in 
recent times. The phenomenon of chaos in low dimensional systems is almost well understood. 
However there is a lot to be investigated in extended systems. Studies on extended systems show 
that the fractal dimension increases linearly with the system size. Further there are measures such 
as multiple positive Lyapunov exponents, spatial correlation length, information functions and so 
on to quantify the spatiotemporal chaos (Tsl [l^ . 

3 Cellular neural/nonlinear networks (CNNs) as reaction- 
diffusion systems 

The most common factor in the autowave process and other pattern formations in reaction- 
diffusion systems is the presence of active nonlinear medium. We have already discussed the 
role of active medium in the previous subsection for the autowave process. However there are 
physical situations where propagation of waves is inhibited beyond certain spatial distances in 
such medium, for example, the failure of electrical impulse in the nerves of patients suffering from 
multiple sclerosis. Unfortunately there are theorems[^ which clearly shows that continuous mod- 
els cannot exhibit such propagation failure. In order to simulate this kind of new phenomenon, in 
addition to the existing several kinds of patterns noted above, it becomes necessary to consider 
discrete versions of diffusively coupled nonlinear dynamical systems to mimic reaction-diffusion 
processes. In many situations, the discrete systems are modeled by appropriate nonlinear elec- 
tronic circuits, for example the cases of impulse propagation along nerve fiber, propagation of 
action potential in cardiac tissues and so on. 

Under these circumstances, it is of great interest to investigate the dynamics of arrays of 
diffusively coupled nonlinear oscillators and systems. In such cases one can often consider an array 
of interconnected locally coupled cells such as neurons, nonlinear circuits, nonlinear oscillators and 
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so on. Such aggregates of cells may be called cellular neural networks (CNNs) in the case of neurons 
and cellular nonlinear networks (again CNNs) in the case of oscillators and circuits jl^, [l9| . 

3.1 Cellular Neural Networks (CNNs) 

In general, a CNN is defined mathematically jl^ by the dynamics of the constituent subsystems 
(state equations of the individual oscillators or cells) and a synaptic law which specifies the interac- 
tion with their neighbours. Such cellular neural networks, namely, interconnections of sufficiently 
large number of simple dynamical units can exhibit extremely complex, synergetic and self or- 
ganizing behaviours. Theoretically one can consider a system of coupled ordinary differential 
equations (odes) to represent a macroscopic system, in which each of the odes corresponds to the 
evolution of the individual subsystem. Such collection or aggregates could represent the discrete 



reaction- diffusion systems\20 . In other words, in reaction-diffusion systems, the CNNs have linear 
synaptic law which approximates the spatial Laplacian operator (nearest neighbour coupling). A 
CNN can be represented by the following four specifications jl^ : 

1. The cell dynamics given by 

Xj = gjiXj) + Ij{xi,X2 ■ . .Xj^i,Xj,Xj + i . . .Xm), j = 1,2 . . .TO, 

where Ij represents to the interaction between the and the remaining cells. 

2. A synaptic law representing the interaction between the cells. For the reaction-diffusion 
CNN this can be given by the discrete spatial Laplacion operator. 

3. Appropriate boundary conditions. 

4. Initial conditions. 

Specific examples are the coupled array of anharmonic oscillators [ pl|, Josephson junctions[p2|, 
continuously stirred tank reactors (CSTR) exhibitingtravelling wavesp3| p4| , propagation of nerve 
impulses (action potential) along the neuronal axon|6|, the propagation of cardiac action potential 
in the cardiac tissues]^ and so on. For the past few years several investigations have been 
carried out to understand the spatiotemporal behaviours of these coupled nonlinear oscillators and 
systems. The studies on these systems include the travelling wave phenomena, Turing patterns, 
spatiotemporal chaos and synchronization!^ |7[ |2|, ^ |3|]. Of particular interest among 

coupled arrays is the study of diffusively coupled driven systems as they represent diverse topics 
like Faraday instability [Q, granular hydrodynamics |l^, self organized criticality[^ and so 
on. Identification of localized structures in these systems has been receiving considerable attention 
very recently |p^ . 

A rather powerful and practical way of studying CNN systems is to model the constituent 
cells in terms of suitable nonlinear electronic circuits, which are then interconnected through 
appropriate linear resistors. The advantages of such arrays of nonlinear electronic circuits is that 
they are quite flexible, that is they can mimic real systems but also can be studied on their own 
merit, they are easy to produce and easy to study experimentally and numerically. From this 
point of view already Perez-Muiiuzuri et al have studied the dynamics of CNNs with the three 
variable Chua's circuit as the basic element and identified several interesting patterns. However 
it will be considerable interest if one investigates the dynamics of very simple diffusively coupled 
driven nonlinear electronic circuits to realize novel spatiotemporal patterns. Recently Murali, 
Lakshmanan and Chua[Q have introduced the simplest second order nolinear nonautonomous 
disspative circuit consisting of a single nonlinear element, namely, the Chua's diode. This simple 
circuit can exhibit a variety of interesting bifurcations, chaos and so on when driven by external 



periodic force ||3^, |37| Therefore it will be of considerable interest to steady the dynamics of 
one and two dimensional arrays of coupled MLC circuts. In this lecture, we wish to give brief 
details of the type of spatiotemporal patterns and other features which arise in arrays of coupled 
MLC circuits with and without the presence of external force. 
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3.2 Murali-Lakshmanan-Chua circuit 



The Murali-Lakshmanan-Chua circuit, Fig. |^(a), is the simplest second order dissipative non- 
autonomous circuit, consisting of Chua's diode as the only nonlinear element . This circuit 
contains a capacitor (C), an inductor (L), a linear resistor (R), an external periodic forcing 
{fsmflt) and a Chua's diode. By applying the Kirchoff's laws to this circuit, the governing 
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Figure 3: (a) Circuit realization of the simple MLC circuit (b) Circuit realization of the Chua's 
diode. The characteristic curve for the Chua's diode {h{x)) for (c) {mo, mi, 7712} — {—1.02, —.55, 
-0.55} and (d) {mo, mi, ma} = {-2.5, 1.5, 0.25} 

equations for the voltage v across the capacitor C and the current iL through the inductor L are 
represented by the following set of two first order nonautonomous differential equations: 

i% = -RiL- RsiL~v + fsinnt, (11) 
at 

where g(v) is a piecewise linear function corresponding to the characteristic of the Chua's diode 
(N) and is given by 

( (' + G2v + {Go-Gi), v>Bp 
g{v) = I e' + Gov, -Bp<v< Bp (12) 

[(' + Giv-{Go~Gi), v<-Bp 

The piecewise nature of the characteristic curve of Chua's diode is obvious from Eq. (|l^). 
The slopes of left, middle and right segments of the characteristic curve are Gi, Go and G2, 
respectively. —Bp and Bp are the break points and e' corresponds to the dc offset in the Chua's 
diode. Rescaling Eq. (|ll|) as w = xBp, iL = GyBp, G = 1/R, u = VlC/G, t = tC/G and e = e'/G 
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and then redefining r as i the foUowing set of normalized equations are obtained: 



X = y-h{x), 

y = —Px — ay + FsiiiLot, (13) 

with 

{e + m2X + {mo - mi) , x > X2 
e + m,QX, xi < X < X2 , (14) 

e + miX — (toq -- mi), x < Xi 

where f3 = (C/LG^), a = [C / LG'^){l + GRs) and F = {ffi/Bp). Obviously h{x) takes the form as 



in Eq. ( |14[ ) with rng Gq/G, mi — Gi/G and m2 = G2/G. The dynamics of Eq. (13) depends 
on the parameters /3, cr, mg, mi, m2, e, w and _F. 

The rescaled parameters in the experimental observations correspond to (3 = 1.0, a — 1.015, 
toq = —1.02, mi = —0.55, m2 — 0.55, e = and lo = 0.75. By varying F one can observe the 
familiar period doubling bifurcations leading to chaos and several periodic windows in the MLC 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

F F 



Figure 4: Bifurcation diagram in the F — x plane (a) for {mo, mi, m2, e, (3, a, lS\ = { — 1.02, 
-.055, -0.55, 0.0, 1.0, 1.015, 0.75} and (b) for {mo, mi, ma, e, /3, a, = {-2.25, 1.5, 0.25, 0.0, 
1.0, 1.0, 0.75} 

circuit. Fig. |4|(a) shows the one parameter bifurcation diagram in the F — x plane \F e (0, 0.7)]. A 
summary of bifurcations that occur in this case for different F values is given in Table |l|. Further it 
is of great interest to consider the parametric choice {mo, mi, ma, e, /?, cr, w} = {—2.25, 1.5, 0.25, 0.0, 
1.0, 1.0, 0.75} which corresponds to the function h{x) having the form as shown in Fig. ^(d). This 
choice of parameters provides the possibility of bistability nature in the asymmetric case in the 
absence of periodic forcing. In this case one can easily observe from numerical simulations that 
the MLC circuit admits only limit cycles for F G (0,0.7). The bifurcation diagram in the F — x 
plane is depicted in Fig. |^(b). 



Table 1: Simrmary of bifurcation phenomena of Eq. (13) 



amplitude (F) 


description of attractor 


0< F < 0.071 


period- 1 limit cycle 


0.071 < F < 0.089 


period-2 limit cycle 


0.089 < F < 0.093 


period-4 limit cycle 


0.093 < F < 0.19 


chaos 


0.19 < F < 0.3425 


period-3 window 


0.3425 < F < 0.499 


chaos 


0.499 < F < 0.625 


period-3 window 


0.625 < F 


period- 1 boundary 
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3.3 Present study 



In the following we discuss the active wave propagation and various spatiotemporal patterns asso- 
ciated with the one and two dimensional CNNs represented by diffusively coupled MLC circuits. 
First we made a critical review on the active wave propagation that occurs in the autonomous 
system by analysing the linear stability properties of the coupled systems. This analysis give us 
a clear understanding of the wave propagation in coupled systems which is essentially the loss of 
stability of the steady states via subcritical bifurcation coupled with the existence of necessary 
basin of attraction for the steady states associated with the coupled system. We also discuss the 
effect of weak coupling on the active wave propagation which causes a blocking of the wave when 
it reaches the weakly coupled cell. Then we look into various spatiotemporal patterns due to the 
presence of Turing instabilty in the absence of external force. 

Further as mentioned earlier it is of great physical interest to study the dynamics of the coupled 
oscillators when individual oscillators are driven by external forces. We study the spatiotemporal 
patterns in the presence of periodic external force and investigate the effect of it on the propagation 
phenomenon and Turing patterns. Depending upon the choice of control parameters, a transition 
from hexagons to regular rhombic structures, hexagons to rolls and then to breathing oscillations 
from hexagons are observed. The presence of external force with sufficient strength removes the 
penta-hepta defect pair originally present in the spontaneously formed hexagonal patterns leading 
to the formation of regular rhombic structures. The inclusion of external periodic force can also 
induce a transition from hexagons to rolls provided there are domains of small roll structures in 
the absence of force. We further show that in the region of Hopf- Turing instability, the inclusion 
of external periodic force with sufficiently small amplitude induces a type of breathing oscillations 
though the system shows a regular hexagonal pattern in the absence of any external force. 

Finally, we also study the spatiotemporal chaotic dynamics of the one dimensional array of 
MLC circuits when individual oscillators oscillate chaotically. In this case, the emergence of 
spatiotemporal patterns depends on the system size. For larger size, above a critical number 
of cells, we observe a controlled space-time regular pattern eventhough the single MLC circuit 
itself oscillates chaotically. However, synchronization occurs for a smaller system size, below the 
threshold limit. 



4 Arrays of Murali-Lakshmanan-Chua (MLC) circuits 

As mentioned in the introduction (see [3.2| ), the circuit proposed by Murali, Lakshmanan and 
Chua (MLC) is one of the simplest second order dissipative nonautonomous circuit having a single 
nonlinear element . Here we will consider one and two dimensional arrays of such MLC circuits, 
where the intercell couplings are effected by linear resistors. 

4.1 One-dimensional array 

Fig. ^ shows a schematic representation of an one dimensional chain of resistively coupled MLC 
circuits. The dynamics of the one dimensional chain can be easily shown to be governed by the 
following system of equations, in terms of suitable rescaled variables, 

Xi = y., - h{xi) + D{xi+i + Xi^i ~2xi), (15) 
Vi = - ay i ~ l3xt + F sin ujt, i = 1, 2, • • ■ , A^, (16) 

where D is the diffusion coefficient, N is the chain length and h{x) is a three segment piecewise 
linear function representing the current voltage characteristic of the Chua's diode and given in 
eq. (|l^). In (|l^ ) mo, mi and m2 are the three slopes. Depending on the choice of mo, mi and m2 
one can fix the characteristic curve of the Chua's diode. Here e corresponds to the dc offset. In our 
analysis, we will consider a few typical forms of h{x) for which the coupld MLC circuits exhibit 
interesting dynamics such as active wave propagation, Turing patterns, spatiotemporal chaos and 
so on. In the following studies we will consider a chain oi N ~ 100 MLC circuits. 
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Figure 5: Circuit diagram showing one dimesnional array of coupled MLC circuits 



4.2 Two-dimensional array 

As in the case of the one dimensional array introduced above, one can also consider a two dimen- 
sional array with each cell in the array being coupled to four of its nearest neighbours with linear 
resistors. The model equation can be now written in dimensionless form as 

= f{xij,yij), (17) 
= -o-yjj - P^ij + D2{yi+i,j + Vi-i.j + Vlj+i + 2/ij-i - iytj) + Fsmujt 
= 9ixz,j,yi,j), (18) 
i,j = 1,2,---,N. 

This two dimensional array has N x N cells arranged in a square lattice. In our numerical study 
we will again take N — 100. 

In the following sections we present some of the interesting dynamical features exhibited by 
the above arrays of coupled MLC circuits such as active wave propagation and its failure, effect of 
weak coupling in the propagation, Turing patterns, effect of external periodic forcing on the Turing 
patterns and spatiotemporal chaotic dynamics. We have used zero flux boundary conditions for 
the study of propagation phenomenon and Turing patterns and periodic boundary conditions for 
the study of spatiotemporal chaos in our analysis. 



5 Spatiotemporal patterns in autonomous CNNs: wave phe- 
nomena and Turing patterns 

Transport processes in living tissues, chemical and physical systems have been often found to be 
associated with a special type of waves, namely, active waves. In the previous section we have 
seen a class of such active waves which arise in many reaction-diffusion systems. The important 
aspect in the reaction-diffusion CNNs is the wave propagation failure which occurs when the 
interconnections are weak, that is, for low values of coupling strength, a feature which can not 
be realized in the continuum limit as proved by Keener in continuous homogeneous reaction- 
diffusion systems. Interestingly, experiments on biological systems shows that there exists a variety 
of situations where such failure in wave propagation do arise. We will identify such possibilities 
in the arrays of coupled MLC circuits. 



5.1 Active wave propagation and its failure in one dimensional CNNs 

To illustrate the wave propagation failure we consider the CNN model described by the array of 
coupld MLC circuits in one dimension, Eqs. (|l^). For this purpose we have numerically integrated 
Eqs. ( p^ ) using fourth order Runge-Kutta method with fixed step. In this analysis we fix the 
parameters at a, mo, wi, m2, e, F} — {1.0, 1.0, —2.25, 1.5, .25, 0, 0} so that the system admits 
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bistability which is a necessary condition to observe a wave front. Zero flux boundary conditions are 
used in the numerical computations, which in this context mean setting xq = xi and x^^i = xm 
at each integration step; similar choice has been made for the variable y also. To start with, we will 
study in this section the autonomous case {F — 0) and extend our studies to the nonautonomous 
case {F ^ 0) in the next section. 

The choice of the values of the parameters guarantees the existence of two stable equilibrium 
points = {cr(mi — mo — e)/(/3 + m2cr), (3{mo — mi — t)/{j3 + m2a)} and P[' = {<7{mi — mo — 
e)/{(i + mia), /3(mo — mi — e)/(/3 + micr)} for each cell, i = 1,2, ... , 100. In the particular case 
corresponding to the above parametric choice, each cell in the array has three equilibrium points 




Figure 6: Space-time plot showing (a) propagation of wave fronts in one dimensional array (100 
cells) of MLC circuits for D = 2.0 and (b) propagation failure for D = 0.4. 

P+ = (3.0,-3.0), Pr- = (-1.5,1.5) and P° = (0,0). Out of these three equilibrium points, P+ 
and P~ are stable while P° is unstable. Due to the asymmetry in the function h{x) for the present 
parametric choices, the basin of attraction of the point P^ is much larger than that of P,[' and it 
is harder to steer a trajectory back into the basin P^ once it is in the basin of P^ . 

Now we choose an initial condition such that the first few cells in the array are excited to the 
P^ state (having a large basin of attraction compared to that of P^) and the rest are set to P[^ 
state. In other words a wave front in the array is initiated by means of the two stable steady states. 
On actual numerical integration of Eqs.(^5|) with N = 100 and with the diffusion coefhcient chosen 
at a higher value, D — 2.0, a motion of the wave front towards right (see Fig. ^(a)) is observed, 
that is a travelling wave front is found. After about 80 time units the wave front reaches the 
100*^ cell so that all the cells are now settled at the more stable state (P^) as demonstrated in 
Fig. ||(a). When the value of D is decreased in steps and the analysis is repeated, the phenomenon 
of travelling wavefronts continues to be present. 

However, below a critical value of the diffusion coefficient {D = Dc) a failure in the propagation 
has been observed, which in the present case turns out to be 13 = Dc = 0.4. This means that 
the initiated wavefront is unable to move as time progresses and Fig. ^(b) shows the propagation 
failure for D = 0.4. 



5.2 Propagation failure mechanism: A case study 

In the above, we have discussed the phenomenon of wavefront propagation and its failure in the 
one dimensional array of MLC circuits by numerically integrating eqs. (|l^) for N ~ 100 cells. In 
order to understand the mechanism behind wave propagation for large coupling strengths and its 
failure for low values of the intcrcell coupling, it is more advantageous to deal with fewer number 
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of cells such as = 3, 4, 5, 6, etc. rather than — 100 because of the diffuculties in handling large 
number of cells and equations analytically. For example one has to deal with the existence of 3^ 
steady states, which is quite large for iV = 100, while it is manageble for A'' = 3, 4 or 5. With this 
fact in mind, we analyse analytically the case iV = 3 and investigate the stability of the stationary 
states in order to understand the nature of bifurcations. A more detailed calculation for A^ = 5 is 
given elsewhere p9|]. One can easily check that for A^ = 3 there are 3'^ = 27 steady states. They 
may be obtained from the following set of defining equations: 

ii = y\~h{x{) ^ D{x2 - x{), 

Vi = -CT%-/3a;j, j = 1,2,3, 

where h{x) is as given in Eq. (|lj) and the parameters are fixed at /3 
mi = 1.5, 7712 = 0.25 and e = 0. 



(19) 

= 1, (7 = 1, mo = -2.25, 



5.2.1 Linear stability analysis 

Now let us consider the wave propagation. It is well known that the propagation failure is due to 
the existence of a large number of stationary states in the system [|l^. Since we have considered 
the coupled system with three oscillators only there exists relatively few steady states. As we are 
looking for the wave front solution, it is sufficient to analyse the stability of only a subset of the 
possible 27 states. This can be done by considering six possible steady states in the following 
way: (i) the ultimate steady states, namely, Xp, Xp and Xp, which exist even in the absence 
of coupling and (ii) the next three steady states which form wave front solutions, namely, X^, 
X^ and X^. The steady states X^, X^ and X^ can be found by assuming the first cell in the 
P+ state and the last cell in the state and allowing the middle cell to be in one of the three 
available states, P+, . For example, X^ state is obtained by considering the middle cell 

in the P+ state. In a similar fashion the rest of the states, X^ and X^ , are obtained. All the six 
steady states can be explicitly given as below: 

X+ = {3,3,3,-3,-3,-3} 
X% = {0,0,0,0,0,0} 

Xp = {-1.5,-1.5,-1.5,1.5,1.5,1.5} (20) 

^ {^Is ' -^28 1 2^3s ) Vis 1 vis J Vzs } 
^% — {.^\si^\s^'^siV\siV%siV%s\ 

= {2;rs>2;2s.a;3^)2/rs>2/2s>%J (21) 



where 



"-Is 



3(81)2 _^ 7Q£, ^25) + _ 3{8D^ + AOD + 25) 
321)2 + 70L» + 25 ' - 32D2 + 701) + 25 ' 



, 1.5(161)2 - 401) - 25) , , 

His ~ ■'^is I ' — -"^i ^: 



X 



32Z)2 + 70Z) + 25 

15(61) - 5) 30Z) 



_ 

161)2 + lOD - 25 ' "^2^ 261)2 ^ loZ) - 25 ' 



^7.5(4Z)-5) 



'Jis ~ -^isJ * — Ij 2, 3, 



" 16Z)2 + 10Z)-25' 

0.6(41)2 - 20Z) - 25) _ ^ 0.3(81)2 ^ 20i:» + 25) 
^ 4L>2 + lOD + 5 ' "^2^ - 4Z)2 + 10L» + 5 ' 

0.3(81)2 _^ 50£) + 25) 
^ 4i:>2 + lOD + 5 ' ^" 



xT^, 7 = 1,2,3. (22) 
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The stability of the above steady states can be found by hnearising (|l^) about them, which leads 
to the following Jacobian matrix: 

( -D-jt, D 10 

D -2D - D Q I 

J= -D--§^ 

-1 0-10 

-1 0-1 

V -10 

where the derivatives, Sx^ '§x~ evaluated at the steady states. The stability is 

guaranteed if all the eigenvalues of the Jacobian matrix (|2^) have negative real parts. One can 
easily find that the states Xp and Xp are asymptotically stable and Xp is unstable irrespective 
of the values of D. However, on a careful analysis on the eigenvalues of the Jacobian for the 
nontrivial steady states, one can find the following: the state Xg is asymptotically stable for 
D < 0.3917 and Xg is asymptotically stable for D < 0.5178 while Xg is unstable irrespective of 
the D values. Thus any initial condition in the neighbourhood of X+ (Xg) will eventually return 
to X+ (Xg) as long as < £> < 0.3916 {0 < D < 0.5177). Further, an arbitrary initial condition 
will evolve to a steady state which lies in the neighbourhood of it. As a consequence, the system 
with a low value of the diffusion coefficient for an initial condition lying in the basin of attraction 
of Xg or Xg will transit to the corresponding stable states X^ or Xg respectively and so a failure 
of propagation occurs. However, for D > Dc, the relevant steady loses stability and consequently 
the systems transits to a suitable final steady state which is stable and this turns out to be Xp. 
As a result wavefront propagation occurs. Further it may so happen that a given initial condition 
may or may not lie in the basin of attraction of a particular steady state Xg or Xg , for a specific 
value of D in which case propagation does occur. Taking this fact also into account, one can prove 
that for the chosen initial conditions described in the numerical analysis above one finally obtains 
D = Dc — 0.4. Thus the wave front propagation occurs when the steay states lose their stability 
via a subcritical bifurcation. 

The above analysis has been extended to the cases = 4, 5 or 6 and one indeed draws the 
same conclusion as above. Thus we can explain the propagation failure at a critical value D = Dc 
as due to the existence of stationary states combined with their basin of attraction. 



\ 



1 




-1 J 



(23) 



5.3 Effect of weak coupling 

In the above, the investigation has been made by considering the system as an ideal one (as far 
as the circuit parameters are concerned). But from a practical point of view, there are defects in 
the coupling parameters which may result in a weak coupling at any of the cell in the array. In 
the following we study the effect of such weak coupling on the propagation of wave front. 

Let us consider a weak coupling at the k^^ cell. By this we mean that the k^^ cell in the array 
is coupled to its nearest neighbour (fc + 1)*^ and (fc — l)^'^ cells by resistors with slightly higher 
values than that of the others. We have studied the effect of this defect on the propagation of 
wave fronts. Numerical simulations have been carried out by considering an one dimensional array 
with 100 cells, where the initial conditions are chosen as in the case of propagation phenomenon in 
regular one dimensional arrays (see Sec. |3.1[ ). From the numerical simulation results, we find that 
there is an abrupt stop in the propagation when the wave front reaches the weakly coupled cell. 
This happens when the coupling coefficient on either side of the k^^ cell in the array has a value 
even above the critical value {D = 0.4) for propagation failure discussed in Sec. 5T. Fig. ^ shows 
a blocking in the propagation of the wave front when the coupling coefficients on either side of the 
k^^ cell {k = 25), which we call as Dfc, is set to 0.47173 with the rest of the coupling coefficients 
set to 1 (D = 1). We observe that the actual blocking occurs when the wave front reaches the fc^^ 
cell. We can say that this is a kind of failure in the propagation because the wavefront will never 
reach the last cell (that is the 100^"^^ cell in the array) by means of blocking. Here the existence 
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Figure 7: Figure showing the effect of weak couphng at the 25 cell in the one dimensional array 
(Eq. (16)). (a) three dimensional space-time plot and (b) density plot. 



of stable steady states at the cell due to the weak coupling is responsible for the blocking of 
wave front. A stability analysis in the present case can be done as in the defect free case and one 
can conclude that the subcritical bifurcation is the cause of blocking. 

5.4 Turing patterns 

Another interesting dynamical phenomenon in the coupled arrays is the formation of Turing pat- 
terns. These patterns are observed in many reaction diffusion systems when a homogeneous steady 
state which is stable due to small spatial perturbations in the absence of diffusion becomes un- 
stable in the presence of diffusion pl|]. To be specific, the Turing patterns can be observed in a 
two variable reaction-diffusion system when one of the variables diffuses faster than the other and 
undergo Turing bifurcation, that is, diffusion driven instability 0. 

Treating the coupled array of MLC circuits as a discrete version of a reaction-diffusion system, 
one can as well observe the Turing patterns in this model also. For this purpose, one has to study 
the linear stability of system (3) near the steady state. In continuous systems, the linear stability 
analysis is necessary to arrive at the conditions for diffusion driven instability. A detailed derivation 
of the general conditions for the diffusion driven instability can be found in Murray]^. For discrete 
cases one can follow the same derivation as in the case of continuous systems by considering 
solutions of the form expi{kj — Ai)[^, Here k and A are considered to be independent of the 
position j (j = 1, 2, • • • , N). For Eq. (3), the criteria for the diffusion driven instability can be 
derived by finding the conditions for which the steady states in the absence of diffusive coupling 
are linearly stable and become unstable when the coupling is present. One can easily show the 
eigenvalues that guarantee the linear stability in the absence of coupling are the roots of the 
characteristic equation 



where fx, fy, Qx and gy are the partial derivatives of / and g in Eq. (3) without coupling coefficients 
{Di — D2 — 0) and evaluated at the steady state. It can be further seen easily that the steady 
state is stable in the absence of coupling if and only if the roots of (4) (Ai and A2) have negative 
real parts. 

Apart from the above condition, in order to satisfy the instability in the presence of coupling 
(Turing instability), atleast one of the roots of the characteristic equation. 



A^ - {fx + gy)>^ + fxgy - fygx = 0, 



(24) 



\l - [k\Di + D2) - {fx + gy)]\s + m{k^) 0, 



(25) 
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with 



m(fc^) = DiD2k'^ - iD2fx + Digy)k'^ + f^gy - fyg^, 

should have positive real part. 

A straightforward calculation shows that the following conditions should be satisfied for the 
general reaction-diffusion system of the form given by Eq. ( p^ : 

fx + 9y < 0, 
fx9v - fy9x > 0, 

fxD2 - gyDi > 0, 

{fxD2 - gyDif - 4DiD2{fxD2 ~ gyDi) > 0. (26) 
The critical wave number for the discrete system (3) can be obtained as 



Combining Eqs. (^6|)-(^7p, one obtains |39|| the condition for the Turing instability such that 

We have applied these conditions to the coupled oscillator system of the present study. For 
this purpose, we have fixed the parameters for the two-dimensional model (3) as {P, a, e,mo, 
Wi,m2} = {0.8, 0.92,0.1, —0.5,0.5,0.5} and verified that this choice satisfies the conditions ( p6| ) 
to (p8|). The numerical simulations have been carried out using an array of size 100 x 100 and 
random initial conditions near the steady states have been chosen for the x and y variables. Figs. 
^(a)-^(d) show how the diffusion driven instability leads to stable hexagonal pattern (Fig. ||(d)) 
after passing through intermediate stages (Figs. ||(a)-||(c)). Further, the spontaneously formed 
patterns are fairly uniform hexagonal patterns having a penta-hepta defect pair. These defects 
are inherent and very stable. 





Figure 8: The spontaneous formation of Turing patterns in an array of 100 x 100 oscillators for 
the parameters f3 = 0.8, <t — 1, mo — —0.5, mi = 0.5, m2 — 0.5, e = 0.1, F = 0.0, uj = 0.75, 
Di^l and D2 ^ 10 in Eq. (3) at various time units (a) T = 50, (b) T = 100, (c) T = 500 and 
(d) T — 2000 (e) the same figure (d) in the Fourier space. 



6 Spatiotemporal patterns in the presence of periodic ex- 
ternal force 

The effect of external fields on a variety of dynamical systems has been studied for a long time 
as driven systems are very common from a practical point of view. These systems can be ei- 
ther spatially modulated or temporally modulated. For example, in a large number of dynamical 
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systems including the DufRng oscillator, van der Pol oscillator and the presently studied MLC 
circuits, temporal forcing leads to a variety of complex dynamical phenomena including bifurca- 
tions and chaos. In particular, it has been shown that the resonant coupling between the forcing 
and the oscillatory modes may lead to several complex dynamical patterns including quasiperi- 
odicity, intermittency and chaos p8| |4l| , |4^ . Also the studies on the effect of external fields in 
spatially extended systems have been receiving considerable interest in recent times ^, p4| . 
Particularly, with the recent advances in identifying localized and oscillating structures and other 
spatiotemporal patterns in driven nonlinear dissipative systems such as granular media, driven 
Ginzburg-Landau equations and so on, it is of special interest to study the effects of forcing on 
arrays of coupled systems such as (3). Motivated by the above, we investigate the effect of external 
forcing on the propagation of wave front and formation of Turing patterns in the coupled MLC 
circuits in one and two dimensions. 



6.1 Effect of external forcing on the propagation of wave fronts 



In this subsection we study the effect of external forcing on the propagation of wave fronts. For this 
purpose we consider an one dimensional array of coupled MLC circuits with initial and boundary 
conditions as discussed in Sec. 5.1. Now we perform the numerical integration by the inclusion 
of external periodic force of frequency uj = 0.75 in each cell of the array (see Eq. (1)). By 
varying the strength, F, of the external force we study the behaviour of t he propagating wave 
front in comparison with the force free case {F = 0) as discussed in Sec. ^.l[ We find that in 
the propagation region {D > Dc), the effect of forcing is just to introduce temporal oscillations 
and the propagation continues without any disturbance (see Fig. ^(a)) as in the case F = (Fig. 
^(a)). Of course this can be expected as the external force is periodic in time. However, interesting 
things happen in the propagation failure region discussed in Sec. 5.1. In this region, beyond a 



certain critical strength of the external forcing, the wavefront tries to move a little distance and 
then stops, leading to a partial propagation. Fig. ^(b) shows such a partial propagation observed 
for F = 0.6 and D = 0.22 (This may be compared with Fig. ||(b)). The phenomenon can be 
explained by considering the propagation failure mechanism discussed in Sec. 5.2 in which one 
may look for a spatially stationary and temporally oscillating wavefront. The initial wavefront 
tries to settle in the nearby stationary state. However, the system will take a little time and space 
to settle due the effect of forcing combined with the transient behaviour of the system. Thus the 
inclusion of external forcing in the propagation failure region can induce the wavefront to achieve 
a partial propagation. 



(a) 



(b) 
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Figure 9: Propagation phenomenon in the presence of forcing: (a) Propagation of wavefront for 
D = 2.0 and F — 0.6 and (b) The partial propagation observed for D — 0.22 and F = 0.6. 
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6.2 Transition from hexagons to rhombs 



It is well known that the defects are inherent in very many natural pattern forming systems. In 
most of the pattern forming systems, the observed patterns are not ideal. For example, the patterns 
are not of perfect rolls or hexagons or rhombs. A commonly observed defect in such systems is 
the so called penta-hepta defect (PHD) pair which is the bound state of two dislocations p5|. 
Experiments on spatially extended systems often show the occurrence of PHD in spontaneously 
developed hexagonal patterns[|6) In the present case also, the existence of PHD pair can be clearly 
seen from Fig. ^(d). In such a situation, it is important to study the effect of external periodic 
force in the coupled arrays of MLC circuits. 




Figure 10: Figure showing the perfect organization of rhombic array in the presence of external 
periodic force with F = 0.25. 

Now we include a periodic force with frequency to and amplitude F in each cell of the array 
and we numerically integrate Eqs.(|l7|) using fourth order Runge-Kutta method with zero flux at 
boundaries. By fixing the frequency of the the external periodic force as w = 0.75 and varying the 
amplitude (F) we analyse the pattern which emerges spontaneously. Interestingly for F = 0.25, 
the defects (PHD pair) which are present in the absence of external force (Fig. ^(d)), gets removed 
resulting in the transition to a regular rhombic array. Fig. ^ shows the gray scale plot of the 
pattern observed for F — 0.25. Thus, from the above we infer that the inclusion of external 
periodic force can cause a transition from hexagonal pattern to rhombic structures. 



6.3 Transition from hexagons to rolls 

In addition to the transition from hexagons to rhombs by the influence of external periodic force, 
there are also other possible effects due to it. To realize them, we consider a different set of 
parametric choice {/3, cr, toq, mi, tos} = {0.734722, 0.734722, -0.874, -0.4715, -0.4715} with e = 
0.15, Di — 1.0 and D2 — 5.0. For this choice the system shows hexagonal patterns with defects 
including domains of small roll structures (Fig. pl](a)). 

Now when the external periodic force is included a transition in the pattern from hexagonal 
structure to rolls starts appearing. By fixing the frequency of the external periodic force again at 
Lu = 0.75, we observed the actual transition from hexagons to rolls as we increase the forcing am- 
plitude (F). Figs. 0(b)-|l^(d) show the gray scale plots for F — 0.15, 0.35 and 0.45, respectively. 
Obviously the transition is due to the existence of small roll structures in the pattern for F — 
which nucleates the formation of rolls in the presence of forcing. 



6.4 Breathing oscillations 

In the above, we have shown that the inclusion of the external periodic force can make a transition 
from one stationary pattern to another stationary pattern like the transition from hexagons to 
rolls. Besides these, are there any time varying patterns? As mentioned above, patterns such 
as localized and breathing oscillations have considerable physical interest. In this regard, we 
considered the parameters a, e, mo, mi, 7712} = {0.734722, 0.734722, 0.10, -0.874, -0.4715, 
—0.4715} with Di — 2.0 and D2 = 5.0 such that a regular hexagonal pattern is observed in the 
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Figure 11: Figure showing the transition from hexagons to rolls for {/?, cr, e, mp, mi, 7712} = 
{0.734722, 0.734722, 0.15, -0.874, -0.4715, -0.4715} with Di ^ 1 and D2 = 5. (a) F = 0, (b) 
F = 0.15, (c) F = 0.35 and (d) F = 0.45. 



absence of external periodic force. From numerical simulations, we observed that a space-time 
periodic oscillatory pattern (breathing motion) sets in for a range of low values of F. Fig. |l^ 
shows the typical snapshots of the oscillating pattern at various instants for the specific choice 
F = 0.05. We have integrated over 10000 time units and the figure corresponds to the region 
T = 4000 — 4014. Typically we find that the breathing pattern repeats itself approximately after 
a period T — 15.0 in the range of our integration. One may conclude that the emergence of 
such breathing oscillations is due to the competition between the Turing and Hopf modes in the 
presence of external periodic force. 



7 Spatiotemporal chaos 

Next we move on to a study of the spatiotemporal chaotic dynamics of the array of coupled MLC 
circuits when individual cells are driven by external periodic force. The motivation is that over a 
large domain of {F, to) values the individual MLC circuits typically exhibit various bifurcations 
and transition to chaotic motion. So one would like to know how the coupled array behaves 
collectively in such a situation, for fixed values of the parameters. For this purpose, we set the 
parameters at {/3, cr, e, mo, mi, m2,cj} = {1.0, 1.015, 0, —1.02, —0.55, —0.55, 0.75}. The uncoupled 
systems exhibit period doubling bifurcations and chaotic dynamics in the presence of external force. 
In our numerical simulations, we have mainly considered the one dimensional array specified by 
Eq. (1) and assumed periodic boundary conditions. 

7.1 Spatiotemporal regular and chaotic motion 

Numerical simulations were performed by considering 50 cells and random initial conditions using 
fourth order Rungc-Kutta method for six choices of F values. The coupling coefficient in Eq. 
( |l5| ) was chosen as D = 1.0. Out of these, the first three lead to period-T, period-2T, period- 
4r oscillations, respectively and the remaining choices correspond to chaotic dynamics of the 
single MLC circuit. Figs. ||(a)-|l|(g) show the space-time plots for F = 0.05, F = 0.08, F = 0.09, 
F — 0.12, F — 0.13 and F — 0.15, respectively. From the Figs. |l3|(a)-|l^(c), it can be observed that 
for F = 0.05, 0.08 and 0.09 the MLC array also exhibits regular periodic behaviour with periods 
T, 2T and 4T respectively in time alone as in the case of the single MLC circuit. However, for 
F = 0.12 and 0.13(Figs. |^(d) and|l|(e)) one obtains space-time periodic oscillations eventhough 
each of the individual uncoupled MLC circuits for the same parameters exhibits chaotic dynamics. 
We may say that a kind of controlling of chaos occurs due to the coupling, though the coupling 
strength is large here. From the above analysis it can be seen that the macroscopic system shows 
regular behaviour in spite of the fact that the microscopic subsystems oscillate chaotically. 
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Figure 12: Snapshots showing breathing oscillations for F — 0.05 and a, e, mo, rni, 7712} = 
{0.734722, 0.734722, 0.15, -0.874, -0.4715, -0.4715} with Di = 2 and D2 = 5 for various time 
units starting from T = 4001. 

Finally for F = 0.15 the coupled system shows spatiotemporal chaotic dynamics (Fig. p^(f)) 
and this was confirmed by calculating the Lyapunov exponents using the algorithm given by Wolf 
et a/[Q. For example, we calculated the Lyapunov exponents for = 50 coupled oscillators and 
we find the largest three exponents have the values Amax — ~ 0.1001, Ai = 0.0776, A2 = 0.0092 
and the rest are negative (see next subsection for further analysis). 



7.2 Size instability, chaos synchronization and suppression of STC 

Since the above study of spatiotemporal chaos involves a large number of coupled chaotic oscilla- 
tors, it is of great interest to analyse the size dependence of the dynamics of these systems. To 
start with, we consider the case of 10 coupled oscillators with periodic boundary conditions and 



numerically solve the system with the other parameters chosen as in Sec. 7.1. The value of F is 
chosen in the range (0.12, 0.15). We find that this set up shows a different behaviour as compared 
to the 50 cells case. Actually the system gets synchronized to a chaotic orbit rather than showing 
periodic behaviour or spatiotemporal chaos as in the case of 50 cells described above. 

To start with we analyse the dynamics for F — 0.12 by slowly increasing the system size 
from = 10. It has been found that the coupled system (^5|) shows synchronized motion for 
iV < 42. This was confirmed by calculating the Lyapunov spectrum which shows only one positive 
exponent with the rest being negative. For example, for N = 42, Amax = 0.1162, with the rest of 
the exponents being negative. The existence of only one positive exponent is a necessary condition 
to have chaos synchronization[|l]. Fig. [I4|(a) shows the dynamics of the 5*^^ cell in the array. 
The system shows entirely different behaviour when we increase the system size to iV = 43. As 



21 



1.5 
1 

0.5 

H 

-0.5 
-1 
-1.5 



(a) 



_l I I L_ 



1.5 
1 

0.5 
- H - 
0.5 
-1 
-1.5 



"T 1 1 r 



(b) 



J I L_ 



1.5 
1 

0.5 
H 
-0.5 
-1 
-1.5 



"T 1 r 



(C) 



J I L 



10 20 30 40 50 10 20 30 40 50 10 20 30 40 50 

i i i 




10 20 30 40 50 

i 



10 20 30 40 50 

i 



10 20 30 40 50 

i 



Figure 13: Space-amplitude plot showing the spatiotemporal periodic and chaotic oscillations in 
50 coupled MLC circuits for various values of external periodic forcing strength: (a) F — 0.05, (b) 
F = 0.08, (c) F = 0.091, (d) F = 0.12, (e) F = 0.13 and (f) F = 0.15. 
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Figure 14: (a) The chaotic attractor at the 5^^ cell for the synchronized state {N = 42) and (b) 
the periodic orbit in the 5^^ cell for the controlled state {N = 43) 



noted in the previous subsection, 7.1, there occurs a kind of suppression of spatiotemporal chaos. 
Fig. |lj(b) shows the resultant periodic orbit in the 5^^ cell of the array. The maximal Lyapunov 
exponent is found to be negative in this case (Amax = —0.001474). Similar phenomenon has been 
observed for F — 0.13 also. 

Next we consider the case of = 0.15 in which the coupled system in the previous subsection 
showed spatiotemporal chaos. From numerical simulations, we again observed a synchronized mo- 
tion for < 31 and the corresponding Lyapunov spectrum shows one positive exponent only with 
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all the other exponents being negative. But for N > 31, the coupled system shows spatiotemporal 
chaos. The Lyapunov spectrum in this case (for N > 31) possesses multiple positive exponents. 
For example, for N = 43, (Amax = Ao = 0.0997, Ai = 0.0633, As = 0.0038). 

The above type of size instability behaviour has also been found in a coupled Rossler system 
as well which undergoes a short wavelength bifurcation |52[. In such cases, one can find the exact 
value of the size below which stable synchronous oscillations occur. However, we find that the 
coupled MLC circuits do not show short wavelength bifurcations. To arrive at a criteria for 
the size instability we computed the Lyapunov dimension. It has been noted that the fractal 
dimension per unit size, dimension density, is an appropriate quantity for the spatiotemporal 



chaotic systems |15|, 16 . A detailed analysis on the size instability and chaos synchronization in 



coupled system (15) will be reported elsewhere. 



8 Conclusions 

In this lecture/article, we have tried to point out how reactive diffusive nonliear systems can give 
rise to a wide variety of spatiotemporal patterns ranging from trivial homogeneous states and trav- 
elling wavefronts to Turing patterns and spatiotemporal chaos. However even wider phenomena 
can be captured by discretized systems in the form of cellular neural/nonlinear networks. As typ- 
ical examples we considered one and two dimensional arrays of coupled MLC circuits and showed 
a wavefront propagation failure can occur due to certain subcritical bifurcations. The onset of 
various interesting patterns including Turing patterns such as hexagons, rhombs, rolls and the 
effect of external forcing on them leading to breathing oscillations were demonstrated. Finally 
how transitions to spatiotemporal chaos and synchronization occurs and the role of system size in 
these transitions have been brought out. 

What has been demonstrated here is only a miniscule of the various phenomena on a space-time 
scale which can occur in nonlinear reactive diffusive systems in general and CNNs in particular. 
Even these basic patterns can play a very useful role in developing synchronized communication 
systems, patterns recognition, image processing and so on. Much more work need to be done in 
order discern all basic patterns which can arise in CNNs and to understand their structure and 
stability. This can in turn will lead to a better understanding of nonlinear dynmaical systems in 
general. 
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